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We apply Density Functional Theory (DFT) and the DFT+U technique to study the adsorption of 
transition metal porphine molecules on atomistically flat Au(lll) surfaces. DFT calculations using 
the Perdew-Burke-Ernzerhof (PBE) exchange correlation functional correctly predict the palladium 
porphine (PdP) low-spin ground state. PdP is found to adsorb preferentially on gold in a flat 
geometry, not in an edgewise geometry, in qualitative agreement with experiments on substituted 
porphyrins. It exhibits no covalent bonding to Au(lll), and the binding energy is a small fraction 
of an eV. The DFT+U technique, parameterized to B3LYP predicted spin state ordering of the Mn 
d-electrons, is found to be crucial for reproducing the correct magnetic moment and geometry of 
the isolated manganese porphine (MnP) molecule. Adsorption of Mn(II)P on Au(lll) substantially 
alters the Mn ion spin state. Its interaction with the gold substrate is stronger and more site-specific 
than PdP. The binding can be partially reversed by applying an electric potential, which leads to 
significant changes in the electronic and magnetic properties of adsorbed MnP, and ~ 0.1 A changes 
in the Mn-nitrogen distances within the porphine macrocycle. We conjecture that this DFT+U 
approach may be a useful general method for modeling first row transition metal ion complexes in 
a condensed-matter setting. 



I. INTRODUCTION 

Metalloporphyrins (of which metal porphines are the 
simplest examples) are stable molecules which exhibit 
a wide range of optoelectronic, magnetic, and me- 
chanical properties^ Deposited and/or self-assembled 
on metal electrodes, porphyrins are attractive candi- 
dates for novel molecular sensors j2 memory^ and light- 
harvesting components'^ Transition metal porphyrins 
are particularly interesting because of the multiple 
spin/electronic states available to them. For exam- 
ple, ligating the metal ions, or simply depositing the 
molecules on electrodes, can preferentially stabilize one 
state versus another. Understanding the detailed molec- 
ular interactions responsible for binding metallopor- 
phyrins to well-characterized metal surfaces is therefore 
a subject of great technological interest, one which has 
received extensive experimental^iSiiiSiSiL&iLiSiAiiiS and 
theoretica l 13 i 14 study. 

Memory and sensor applications of transition metal 
porphyrins rely on electrochemically-induced switching 
of the electronic and magnetic statesiS^ The field depen- 
dent self-assembly and mobility of porphyrin molecules 
on electrodes have also been the subject of experimen- 
tal interest under ultra-high vacuum conditions^ and 
in solution*^ Electric-field induced changes in the elec- 
tronic states can lead to structural changes in porphyrin 
molecules, and this may impact their functions and 
reactivity in subtle ways^ Macrocycle conformational 
changes in nickel and zinc porphyrins have recently been 



the subject of a series of investigations>i§iiLi£ Attach- 
ing ligands to Ni(II) or photoexcitation can switch the 
low-spin nickel ground state with a (d z 2) 2 configuration 
to the high-spin expanded {d z 2){d x 2_ y 2) state, so that 
core expansion causes the porphyrin ring to undergo a 
conformational change from a ruffled to a dome-like ge- 
ometry. When outlying porphyrin appendages or bri- 
dles are attached to the Ni porphyrins, this small but 
energetic motion triggered by nickel spin state changes 
can potentially be harvested as nanomechanical motion 
in molecular switchesiiSiii It would be of great interest 
to extend this general principle to other transition metal 
porphyrins and related molecules deposited on gold elec- 
trodes. In such cases, the transition metal charge or spin 
states can be altered by electrochemical means, and the 
resulting conformational changes can be monitored using 
atomic force microscopy. 

Given the wide range of interest in porphyrins ad- 
sorbed on metal surfaces, it is important to elucidate 
the precise nature of metalloporphyrin adsorption on and 
interaction with metal electrodes. Experimentally, it is 
known that certain substituted porphyrins and related 
molecules lie flat on Au(lll) surfaces p^iSaLiSiii Even 
when they contain bulky substitutents, the molecules 
are only slightly distorted from a planar adsorption 
geometry^ It is also possible to have substituted por- 
phyrins self-assembled into a tube-like geometry and lie 
edgewise on material surfaces if they are tethered using 
sulfide or other linkages 4£ 

Density functional theory (DFT) might appear to be 



the theoretical method of choice to shed light on the 
binding energies, spin states, geometries, and external 
electric field effects of porphyrins adsorbed on metal 
surfaces. DFT simultaneously addresses the electronic 
and geometric properties of the composite molecule- 
metal system. It is a formally exact method, but 
the quality of its predictions depends on the approx- 
imate exchange correlation functional used. An early 
DFT calculation 1 ^ was performed prior to detailed ex- 
perimental studies of porphyrin adsorption on atom- 
istically flat gold surfacesi WV 2 (For the purpose of 
this work, we ignore the herring-bone surface recon- 
struction on Au(lll)ji£ which occurs on long length- 
scales.) Motivated by molecular electronics experiments 
on self-assembled molecules which are somewhat simi- 
lar to porphyrins^ Lamoen et al^- examined palladium 
porphines (PdP) adsorbed edgewise on Au(lll), and re- 
ported a large, ~ 10 eV, binding energy. This calcula- 
tion was performed using the local density approximation 
(LDA), which tends to overestimate the binding energy. 
A subsequent DFT worki^ considered PdP adsorbed in 
a similar edgewise geometry on Al(lll) surfaces. This 
study reported a sub eV adsorption energy using LDA, 
and an even smaller binding energy using the Perdew- 
Burke-Ernzcrhof (PBE) 21 version of the generalized gra- 
dient approximation (GGA) to the exchange correlation 
functional. GGA is generally more accurate than LDA 
for treating molecule-material binding 2 ^ and has been 
widely used in such studies^ 

Guided by recent experiments and building on these 
earlier DFT results, one of the goals of this work is to 
apply the PBE functional to re-examine PdP adsorption 
on the Au(lll) surface. We consider both the flat adsorp- 
tion geometry found in experiments where porphyrins are 
not tethered to metal surfaces (SiiSiiiiiS and the edgewise 
geometry adopted in earlier calculations. 

PdP can be considered a prototypical metallopor- 
phyrin which adsorbs weakly on Au(lll). For many 
transition metal compounds and materials, such as sec- 
ond row elements like Pd, most implementations of DFT 
appear to yield the correct ground spin state and elec- 
tronic properties, and they are valuable techniques to 
apply. 2 ^ However, we are also interested in first row tran- 
sition metal porphyrins (e.g., Ni(II) and Mn(II), which 
are pertinent for studying effective ion size-driven con- 
formational changes in porphyrin molecules) ^fiiiLiS Here 
we encounter a dilemma that has hindered DFT studies 
of composite systems involving first row transition metal 
ions. In brief, DFT implementations with non- hybrid 
exchange correlation functionals such as PBE 2 — and 
PW91, 24 widely used in the condensed matter physics 
community, treat gold and other metal surfaces accu- 
rately but are less accurate for interactions between lig- 
and/crystal field and first row transition metal ions. The 
Becke 3-parameter Lee- Yang-Parr (B3LYP) 2 ^ hybrid ex- 
change correlation functional correctly predicts the high- 



spin ground state for MnP, but is at present too costly 
for modeling metal surfaces 

Accurate predictions for bulk transition metals lat- 
tice constants and magnetic properties were among the 
early successes of non-hybrid GGA, 29 PW91 and PBE 
are among this class of GGA. However, these function- 
als underestimate the exchange interaction among the 
strongly localized, partially filled 3d orbitals in first row 
transition metal ions. This leads to underestimation of 
the stability of high-spin states in some first row transi- 
tion metal ions in impurity centers or ligand fields, where 
non-hybrid GGA often predicts ground states with in- 
correct spin multiplicitiesi22*2ii2 2 i Indeed, a recent ab ini- 
tio molecular dynamics simulation of Mn(II) centers was 
compelled to apply a high-spin constraint throughout the 
trajectory^ 

Hybrid exchange correlation functionals contain frac- 
tional non-local or Hartree-Fock exchange, which 
amounts to 20% for BSLYPp^ This admixture is appar- 
ently the right amount to reproduce the experimental 
spin ordering and energy splitting between spin states in 
many complexes between ligands and Mn^SS*^ as well 
as other first row transition metal ions^& and is widely 
and successfully used for Mn centersi 34 i 3 ^' 3 9i 3 7i 38 i 3 ? A 
slightly smaller admixture yields better results in other 
cases i 3 ^ 3 ^ At present, applying these hybrid functionals 
in simulation cells with periodic boundary condition is 
computationally costly — up to 100 times more so than for 
non-hybrid functionals^ This factor may further depend 
on system size. As such, B3LYP has only seen prelimi- 
nary applications in condensed phase systems 4& To our 
knowledge, it has not been applied to the large, slab-like, 
periodically replicated supercells needed to study por- 
phyrin adsorption on metal surfaces^ or for that mat- 
ter, for large scale DFT calculations that require periodic 
boundary conditions, such as in aqueous systems via ab 
initio molecular dynamics^ 3 - As will be shown, periodi- 
cally replicated simulation cells with adequate Brillouin 
zone sampling are crucial for modeling porphyrin adsorp- 
tion on gold surfaces. 

We note that hybrid functionals are not universally 
successful for all transition metal species. For exam- 
ple, (1) hybrid functionals tend to overestimate the sta- 
bility of TVs 1 (TV — l)d n+1 electronic configurations over 
Ns 2 (N — l)d n ones in transition metal atoms^ (So does 
non-hybrid GGAml) This consideration does not play a 
role when the metal s-orbital is no longer available due 
to its involvement in doubly occupied bonding orbitals, 38 
or in Mn(II)P and Mn(III)P, which have largely empty 
Mn 4s orbitals. (2) The binding energies and bond 
lengths of many transition metal dimers are more ac- 
curately predicted using non-hybrid rather than hybrid 
GGAi^ This issue may be related to (1). Thus, com- 
plexes with Mn-Mn bonds may be particularly challeng- 
ing for hybrid functionals 4244 (3) For as yet unknown 
reasons, B3LYP is ambiguous regarding the stability of 
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the high-spin electronic configurations of some first row 
transition metal porphines, such as ligated FeP^* 4 ^ while 
non-hybrid GGA performs considerably worsej^ 5 - a larger 
admixture of Hartree-Fock exchange seems necessary in 
this case. 

A promising alternative method to characterize the 
correct spin state is to use Quantum Monte Carlo (QMC) 
methods, which up to now have been used to study 
porphine (H2P) which lacks coordinated metal ions4£ 
Another option is to impose the experimentally known 
spin polarization in the entire supercell^ As will be dis- 
cussed, however, this approach may not yield accurate 
metal-nitrogen distances within the porphyrin ring. Fur- 
thermore, when a metal electrode is present, the bulk 
metal valence electrons can be excited to the conduction 
band to give arbitary spin multiplicities at negligible en- 
ergy penalties. So, a global spin constraint will likely 
change the bulk metal magnetic behavior without affect- 
ing the adsorbed transition metal ion complex. Clearly, 
new computational techniques are needed to treat por- 
phyrin molecules adsorbed on metals. 

In this work, we apply the DFT+U method 48 ' 49 to 
treat the composite manganese porphine-gold metal sys- 
tem. This method, which emphasizes the role of on-site 
screened coulomb interactions, has successfully predicted 
the correct electronic ground state for transition metal 
oxide crystals such as NiO£S and L&CoO^r^ where tra- 
ditional non-hybrid GGA has failed. This mean field ap- 
proach augments the DFT exchange correlation interac- 
tons among electrons in a subset of orbitals (Mn 3d in our 
case) with Hartree-Fock-like interactions parameterized 
with coulomb (U) and exchange (J) terms. Electrons in 
these orbitals still interact with the rest of the system 
(including the C, H, N atoms in the porphine ring and 
the Au(lll) substrate) via the DFT formalism. The spe- 
cific implementation we use is described in more detail 
in the Supporting Information section. 

Formally speaking, the parameters U and J are related 
to coulomb and exchange integrals pSLS In practice, U 
has been set at values required to achieve agreement with 
experiments. 53 When dealing with molecules embedded 
with transition metal ions, we propose that U and J can 
be fitted to either experiments or predictions via other 
gas phase theoretical methods such as high level quantum 
chemistry, quantum Monte Carlo, and hybrid functional 
DFT — as long as the latter are known to agree with ex- 
periments. The DFT+U technique is readily amenable to 
periodically replicated simulation cells, and is thus well 
suited for studying metal porphine molecules adsorbed 
on gold surfaces. We note that this method is empiri- 
cal in nature, and it assumes that correlation between 
the 3d electron can be approximated as static* 5 ^ Never- 
theless, it has seen wide and successful applications to 
many solid state materials where LDA/GGA treatments 
fail for electronic/magnetic properties ^^0^^ ft may 
also become a general technique for ab initio molecular 



dynamics study of first row transition metal centers in 
condensed phases 

As a proof of principle, we will focus on a second pro- 
totype porphyrin: manganese (II) porphine (MnP). MnP 
exhibits behavior drastically different from that of PdP 
discussed above. PdP is umambigously stable in the low 
spin state, and PBE and DFT+U methods yield simi- 
lar predictions. However, Mn(II) complexes in general 
and Mn(II) porphyrins in particular are known to have 
high-spin gr ound states. 56 ! 57 ' 58 We will show that PBE 
incorrectly predicts the intermediate spin state to be the 
MnP ground state. On the other hand, by fitting U to 
B3LYP results (which yield the experimental high spin 
ground state), both the high/intermediate spin energy 
splitting and the gas phase MnP geometry are success- 
fully reproduced using the DFT+U method. 59 

Having validated the DFT+U technique for isolated 
MnP, we further apply it to investigate the geometric, 
electronic, and magnetic properties of MnP adsorbed on 
Au (111) surfaces. We also examine the effect of an elec- 
tric field on the adsorbed MnP, and show that it can 
induce significant changes in these properties. The MnP 
structure is particularly strongly affected by the applied 
field, with the Mn-N distance in the porphine ring in- 
creasing by 0.13 A. This suggests that Mn porphyrins 
may exhibit conformational changes related to those seen 
for Ni porphyrins. 16 

To summarize, we have conducted a comprehensive 
study of the adsorption of two prototype transition metal 
porphine molecules, PdP and MnP, on gold surfaces. The 
effect of an applied electric field on the properties of ad- 
sorbed MnP is investigated. We also emphasize the suc- 
cess and importance of applying the DFT+U technique 
to treat the first row transition metal ion Mn(II), and 
conjecture that this technique may be widely applicable 
to condensed phase systems. 

This paper is organized as follows. Section II describes 
the methods and models used in our calculations. Sec- 
tion III details the results on isolated porphine molecules 
and porphines adsorbed on Au(lll). Section IV con- 
cludes the paper with further discussion of the results 
obtained. 



II. METHODS AND MODELS 

DFT /PBE and DFT+U calculations are performed us- 
ing the Vienna ab initio Simulation Package (VASP), 60 
the projector augmented waves method^ and associ- 
ated pseudopotentialspS For details of the DFT+U im- 
plementation, see Bengone et alr& and the Supporting 
Information. The cutoff for wave functions is set to 
400 eV. LDA calculations are performed with both VASP 
and SeqQuest) 6 ^ using ultrasoft pseudopotentials^ 5 - and 
norm conserving pseudopotentials, respectively. An en- 
ergy convergence criterion of 10 -3 eV for each atomic 
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FIG. 1: A metal-porphine molecule, (a) without ligands, 
(b) attached to a ligand (Cl _ in this side view example). 
Carbon, nitrogen, hydrogen, metal, and chlorine atoms are 
depicted as grey, blue, white, red, and green circles, respec- 
tively. The metal-nitrogen (dn), metal-porphine plane (d p ), 
and metal-ligand (d\) distances are illustrated in the figure. 
As depicted, the metal ion in panel (b) is displaced towards 
the ligand atom, and d p is defined as positive ("above the 
macrocycle ring"). The Au(lll) substrate, if present, is be- 
low the macrocycle on the oposite side to the ligand. 



are found with Ahlrichs and co-workers' all-electron TZV 
basis^ 5 . on the metal ion and the 6-311G** basis for all 
other atoms, followed by a single point energy calcula- 
tion using the Wachters-HaySi all electron basis set aug- 
mented with a diffuse d-function, 6-311+G, on the metal 
ion, and 6-31G* on all other atoms. 

Molecular mechanics calculations are performed with 
the LAMMPS molecular dynamics software^ modi- 
fied to utilize the porphyrin potential of Shclnutt and 
coworkers— The gold substrate is 6 layers thick, and has 
dimensions 41 A in the x and y directions. Au atoms are 
held fixed throughout the simulations; the gold-porphinc 
interactions are calculated using previously published 
force fields— Molecular mechanics total energy mini- 
mization calculations are converged to sub-meV levels. 



III. RESULTS 



configuration is enforced. 

Calculations on isolated porphyrin molecules (Fig. 
are performed using 16x16x10 A 3 simulation cells. For 
porphines with metal ions ligated with H2O or Cl~ 
groups, cell sizes up to 16 x 16 x 15 A 3 are used, and dipole 
corrections are applied to remove the coupling between 
periodic imagesi^SiSLS These simulation cells converge 
the spin splittings to 10 meV or better. 

To model porphines deposited flat on Au(lll) surfaces, 
we use supercells of lateral size 14.79x15.37 A 2 . They 
contain three to six gold layers of 30 Au atoms each, and 
we apply various Brillouin sampling schemes. The cell 
lengths in the z-direction are 20, 22, 25, and 28 A when 3, 
4, 5, and 6 layers of gold Au are present, respectively. The 
Au atoms in the bottom layer are fixed at their bulk face 
center cubic (FCC) positions— *22i Dipole corrections are 
also applied for these slab geometry calculations 

To study edgewise adsorption of PdP on Au(lll) us- 
ing LDA, we adopt the supercell of Ref. 13, except that 
the cell dimension in the z-direction is increased to 27 A 
to minimize the coupling between periodically replicated 
images. The supercell contains three gold layers totaling 
36 Au atoms, and a PdP molecule with two edge protons 
straddling a gold atom (see Fig. [2>). The equilibrium 
PBE lattice constant is larger than that of LDA, and 
PBE investigations of this adsorption geometry employ 
a simulation cell of size 5.91x15.36x27 A 3 . Monkhorst- 
Pack grids^i of density up to 4x2x1 are used to sample 
the Brillouin zone. 

Optimizations and single-point energy calculations 
of gas phase porphinc molecules are computed with 
the program Jaguar— using the DFT B3LYP hybrid 
functional 2 ^ with the 6-311G** basis set for second row 
atoms, and the lacv3p** basis set for Mn that includes 
Effective Core Potentials (ECPs). 73 Additional optimiza- 
tion calculations are performed using all-electron basis 
sets with the Gaussian 03 code— Optimized energies 



A. Isolated porphyrin molecules 

This section describes the geometries of isolated por- 
phine molecules, and the energy splitting as a function 
of the spin multiplicity, as predicted by the different 
exchange correlation functionals. All symmetries are 
turned off in PBE and DFT+U calculations, which are 
performed using the VASP code and are converged to 
a few meV. B3LYP calculations for MnP enforce D2h or 
D4h symmetries, depending on the Jahn- Teller distortion 
dictated by the degeneracy. Upon binding to a ligand or 
in the presence of a substantial electric field, the metal 
ion displaces out of the porphine ring towards the ligand 
atom by a distance d p , where the porphine ring position 
is taken to be the average z coordinate of the 36 porphinc 
atoms excluding the metal ion. The distance between the 
transition metal ion and the ligand atom closest to it is 
defined as d\. These definitions are illustrated in Fig.^ 

Table ^ shows that B3LYP predicts a high spin (S = 
5/2) ground state for MnP, more stable by 0.21 (0.25) eV 
over the intermediate spin-state (S = 3/2). The val- 
ues without/ with parentheses are computed using the all 
electron (TZV/6-311G*7/6-311+G/6-31G*) basis set, 
or lacv3p** (Mn) and 6-311G** (other atoms), respec- 
tively. They show that neither the high spin-intermediate 
spin energy splitting, nor the molecular geometry, is 
strongly affected by the basis set used. (We have listed 
the average of the two Mn-N distances in the case of 
quartet MnP, which has a D2h symmetry.) The stabil- 
ity of the high spin state is in agreement with experi- 
ments on substituted manganese porphyrins. 5 ^ PBE, on 
the other hand, incorrectly predicts that S — 3/2 is 
more stable than S = 5/2 by ~ 0.50 eV. When we set 
U = 4.2 eV, J = 1.0 eV within the DFT+U approach, 
the 5 = 5/2 spin polarized state becomes more stable 
than the S = 3/2 state by 0.23 eV. This is sufficiently 
close to the B3LYP results for our purpose— — Since 
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TABLE 1: Spin states and relative energies of metal por- 
phines, and some intramolecular distances: metal-nitrogen 
(<1n), metal-porphine ring (d p ), and metal-ligand (di). Only 
the ground state of Mn(H20)P is tabulated. Energies and 
distances are in units of eV and A, respectively. B3LYP re- 
sults with/ without parantheses are computed using different 
basis sets (see text). 
Experimental results in the solid phase: 
"tetraphenylporphinatopalladium (II) (PdTPP); 85 
b octaethylporphinatopalladium (II) (PdEOP)^ 
c tetraphenylporphinatomanganese (II) (MnTPP); 57 
d tetraphenylporphinatomanganese (III) 
chloride (MnCTPP)i^ 

DFT+U predicts that all Mn(II)P and Mn(III)P species 
examined in this work are high spin, these ground state 
predictions are not affected by spin contamination^ We 
emphasize that we use B3LYP as a standard to param- 
eterize our DFT+U work because it yields the correct 
high-spin ground state for MnP(££ as discussed in the 
introduction, B3LYP is not universally successful for all 
transition metal compounds. 

With these parameters, the ground state Mn-N dis- 
tance is predicted to be cZn = 2.090 A, very similar to 
the B3LYP prediction, and within 0.005 A of the experi- 
mental value of 2.085 A, measured for the corresponding 
tetraphcnylporphyrin derivative^ On the other hand, 
even when the total spin is constrained to S — 5/2, PBE 
predicts a Mn-N distance of 2.051 A, or 0.034 A smaller 
than the experimental value. Thus, even with a high- 
spin constraint, PBE underestimates Mn-N distances in 



Mn(II)P, while the DFT+U approach can be parameter- 
ized to yield both the correct spin state and a reasonable 
metal-porphine geometry. 

Attaching a H2O ligand to Mn does not qualitatively 
change the above analysis, with PBE still predicting an 
incorrect S = 3/2 spin state and a Mn-N distance sub- 
stantially smaller than that predicted using the DFT+U 
method. 

The energy splitting between the high and intermedi- 
ate spin states is sensitive to U, while the Mn-N distance 
is less sensitive. For example, setting U = 3.8 eV instead 
of 4.2 eV yields a 0.12 eV splitting, but the Mn-N bond 
length remains essentially unchanged at 2.088 A. 

In MnPCl, spectroscopic measurements have con- 
firmed that Mn is in the high spin Mn(III) oxidation 
stated in agreement with DFT+U and PBE predictions. 
These methods also predict similar Mn-N bond lengths 
that overestimate the experimental value in substituted 
porphyrins by up to 0.04 ASS4 This discrepancy may 
be related to the ~ 0.1 A difference in the predicted 
and measured out-of-plane Mn displacement. In the gas 
phase, the predicted change in Mn-N distance between 
the Mn(III) and Mn(II) oxidation states is thus underes- 
timated. Despite this, as will be seen, DFT+U predicts 
a significant increase in this bond length when Mn(III)P 
adsorbed on Au(lll) is switched to a Mn(II)-like oxida- 
tion state with an electric field. 

PBE predicts that PdP is low spin, S = 0, in 
agreement with experiments^^ So does DFT+U (not 
shown). In general, the 4<i electrons in second row transi- 
tion metal ions are less localized than 3d electrons in first 
row transition metal ions, and non-hybrid GGA methods 
such as PBE appear more reliable in predicting the order- 
ing of their spin states than is the case with 3d electron 
systems^ 



B. Normal-Coordinate Structural Decomposition 

The conformations of the porphine macrocycles in the 
predicted structures are further examined using Normal- 
Coordinate Structural Decomposition (NSD). 89 NSD has 
emerged as a useful method for analyzing deformations 
of tetrapyrrole macrocycles in heme proteins and in syn- 
thetic proteinsj2S and such analyses are pertinent to the 
potential use of transition metal porphyrins adsorbed on 
metal electrodes as conformational switchesiiSiii Com- 
plete NSD analyses of selected structures from Table I 
are given in Tables S4-S13 of the Supporting Information. 
The PBE predicted structure for PdP and DFT+U pre- 
dicted ones for MnP and MnClP show negligible amounts 
(less than 0.03A) of nonplanar deformation, consistent 
with the crystallographic data for other transition metal 
porphines (see Table S13). For comparison, highly sub- 
stituted and very nonplanar porphyrin macrocycles ex- 
hibit up to 4 A deformations in the soft ruffling (Bi u ) or 
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saddling (B2 U ) modest 

NSD also provides details of the in-plane deformations 
present in the isolated porphine macrocycles^ In-plane 
deformations are not usually analyzed in any detail be- 
cause they are naturally intertwined with out-of-plane 
deformations in substituted porphyrins^ making it dif- 
ficult to isolate the former. However, in the case of the 
porphines where the macrocycles are all nominally pla- 
nar, it is possible to see the in-plane deformations re- 
sulting from the changes in the size of the metal com- 
plexed to porphine. For example, in the structures of 
PdP (PBE), MnClP (DFT+U), and MnP (DFT+U), the 
metal-nitrogen distances are 2.035, 2.045, and 2.090 A, 
respectively. (The X-ray structure of nickel (II) por- 
phine exhibits an even smaller metal-nitrogen distance 
of 1.951 A because of the small size of the Ni(II) ion; 
see Table S13.) The NSD studies reveal a correspond- 
ing increase in the first-order Ai g (breathing) deforma- 
tion as the porphine ring expands to accommodate the 
larger metal ion. The first order Ai g deformations are 
+0.15, +0.11, and +0.37 A (and -0.16 A for NiP). Fi- 
nally, we note that the discrepancy between the metal- 
nitrogen distances predicted for PdP and measured in 
PdTPP and PdOEP may be related to the in-plane defor- 
mations that accompany the out-of-plane deformations 
in the latter, as the crystal structures of these Pd com- 
plexes adopt nonplanar conformations of the type that 
are known to shorten the metal-nitrogen bondZ& (e.g., the 
large ruffling (Bi„) out-of-plane deformation in PdTPP 
(see Table S5)). 

C. Palladium (II) porphine on Au(lll) 

We consider two adsorption geometries: a flat adsorp- 
tion geometry suggested by experiments on substituted 
porphyrins&iSiiiii^ (Fig|2Ji); and the edgewise porphyrin 
stack configuration examined in Ref. 13, with successive 
macrocycles separated by ~ 5.9 A (Fig[2>). 

1. Edgewise adsorption 

When we attempt to adsorb PdP on Au(lll) in the 
edgewise geometry depicted in Fig. |5Jd using the PBE 
exchange correlation functional, we find little or no at- 
traction between the PdP tube and the gold substrate. 
In fact, PBE even predicts that the PdP molecules re- 
pel each other when they reside 5.91 A apart in the self- 
assembled tube-like geometry in the absence of gold. The 
formation energy of this PdP tube is +0.06 eV when us- 
ing 4x2x1 Brillouin zone sampling. The small repulsion 
is partly due to the fact that PBE underestimates van 
der Waals attractions. Surprisingly, we find that the Y- 
point Brillouin zone sampling used in Refs. 13 and 14 
yields an (unconverged) 0.38 eV repulsion. This is de- 
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FIG. 2: Palladium porphine on Au(lll). (a) Flat adsorp- 
tion geometry, with Pd on the top site. The top, hollow, and 
bridge sites are illustrated, (b) Edgewise adsorption. Striped 
circle: N; black: Pd; dark grey: C; white: H; light grey: 
Au. Periodic boundary conditions are used in all three direc- 
tions (a vacuum layer separates the gold surfaces); only the 
gold atoms in the primitive simulation cell are shown. For 
the geometry in (b), the porphine molecules are in fact self- 
assembled in a tube-like stack along the gold surface. 



spite the fact that there is at most a 0.04 eV dispersion 
in the PdP valence electron states near the Fermi level. 
Regardless of which PdP reference energy is used — the 
isolated porphine or porphine tube — we find no attrac- 
tive interaction between PdP and Au(lll), in contrast 
to the ~ 10 eV binding energy as previously reported by 
Lamoen et ali&. 

Lamoen et al. applied LDA, not the PBE exchange 
correlation functional. They reported that the closest 
distance between PdP edge protons and Au atom is 
1.78 A, suggesting a strong covalent bond between H 
and Au. The large adsorption energy was rationalized 
by comparing this system with the interaction of hy- 
drogen molecules on transition metal surfaces, 13 While 
LDA is known to predict over-binding of molecules on 
metal surfaces, the discrepancy between our predictions 
and Ref. 13 is too large to be explained by the differ- 
ent treatment of exchange/correlation. To resolve this 
discrepancy, we revisit the edgewise adsorption binding 
geometry using LDA. We apply both the VASP code^S 
with ultrasoft pseudopotentials and a plane wave basis, 
and the SeqQuest code^i which applies norm-conserving 
pseudopotentials and localized Gaussian basis sets. We 
use two different DFT packages to ensure that our predic- 
tions are not artifacts of the pseudopotentials employed. 
With both of these DFT codes, our predicted equilibrium 
Au lattice constant of 4.078 A agrees with Ref. 13. How- 
ever, starting with a geometry where two edge PdP pro- 
tons are 1.8 A away from surface Au atoms, PdP expe- 
riences large forces and relaxes away from the substrate. 
The optimized distance of closest approach between PdP 
protons and Au atoms is 2.36 A, and is associated with 
an adsorption energy of only 0.25 eV per PdP molecule. 
(The relaxed geometries are listed in Supporting Infor- 
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mation section.) Sampling different adsorption sites by 
moving the PdP center of mass in the x-y plane results 
in only small (~ 0.03 eV) variations in -Ebind- 

We conclude, based on the results of two very different 
DFT codes, that the attraction between Au(lll) and 
edgewise adsorbed PdP is weak, of the order of 0.25 eV 
within the LDA approximation, not the ~ 10 eV reported 
previously^ Note that a DFT work which adopts the 
edgewise adsorption geometry of Ref. 13 to study PdP on 
Al(lll) also reported a LDA binding energy of a small 
fraction of an eV>i^ 

We emphasize that this 0.25 eV attractive interaction 
is computed using LDA. In the remainder of this work, 
we will use PBE or DFT+U based on PBE. The PBE ex- 
change correlation functional predicts a negligible bind- 
ing energy between an edgewise adsorbed PdP tube on 
Au(lll). 



2. Flat adsorption geometry 

Next, we examine a flat adsorption geometry. We 
place the center of mass of the PdP molecule, namely 
the Pd atom, atop the hollow, bridge, and top sites of 
the Au(lll) surface, and carry out geometry optimiza- 
tion. L-point sampling and the 3-layer, 90 atom Au(lll) 
supercell described previously are used. Due to the weak 
interaction between PdP and Au(lll), geometric relax- 
ation along the ^-direction is slow. Nevertheless, when 
PdP previously optimized in the gas phase is initially 
placed 3.17 A from the top layer of Au(lll) atoms, the 
porphine ring eventually relaxes to a distance ~ 3.5 A 
above the gold surface. Here, the z coordinate of the top 
layer of the gold substrate is averaged over all atoms in 
that layer, and the porphine ring z-position is averaged 
over all porphine atoms other than Pd. 

The binding energies for all three sites are similar, and 
are between 0.255 eV and 0.270 eV (Table They are 
consistent with a weak, non site-specific van der Waals 
interaction between PdP and Au(lll). There are mini- 
mal changes in the PdP geometry compared to the gas 
phase system, with the Pd-N distance expanding by less 
than 0.01 A, and the Pd(II) ion displacing towards the 
gold substrate by less than 0.1 A. PBE underestimates 
van der Waals forces between PdP and Au. For a quali- 
tative comparison, our molecular force field calculation 78 
yields a 0.54 eV binding energy for bare porphine (H2P) 
molecules adsorbed in a flat geometry on Au(lll). In- 
creasing the number of gold layers from 3 to 4 increases 
the PBE i^bind by only 0.005 eV, showing that the calcu- 
lation is adequately converged with respect to the system 
size employed. 

The absence of PdP-Au(lll) covalent bonding is con- 
firmed by examining the electron density distributions 
and the electronic density of states (DOS). The electron 
densities associated with PdP and the Au surface pre- 



Site £bind d p dAu rf p /Au 

top -0.270 2.040 -0.079 3.40 3.48 
bridge -0.267 2.039 -0.079 3.40 3.48 
hollow -0.255 2.039 -0.064 3.39 3.45 

TABLE 2: Flat adsorption geometry PdP binding energies 
on various Au(lll) surface sites and some intramolecular dis- 
tances: metal-nitrogen (dn), metal-porphine ring (dp), metal- 
gold surface (dAu), and porphine- gold surface (d p /A u )- Calcu- 
lations are performed using 3-layers of Au atoms and T-point 
sampling. Energies and distances are in units of eV and A, 
respectively. 

sented in Fig. [3^ show minimal overlap. Figure \S]p su- 
perimposes the DOS of an isolated PdP (inverted curve 
in the upper half of the figure) with that of the adsorbed 
PdP-Au(lll) complex (lower curve). We align the two 
sets of DOS by performing calculations where PdP is 
moved far from Au(lll) surface, and line up the lowest 
occupied orbitals, namely the Pd Ap states. We find that 
a 6 or 8 A separation is sufficient to converge the DOS 
alignment to the infinite separation limit. The DOS's 
show that the Fermi level of the PdP-Au complex resides 
within the highest /lowest occupied molecular orbital gap 
of PdP, far from PdP d-orbital levels. Thus, we do not 
expect any charge transfer between PdP and Au. Fig- 
ure Efb) also depicts the orbitals in isolated PdP and 
the PdP-Au complex which exhibit substantial Pd Ad 
character. There are more than 5 such states because 
of hybridization between the Pd and nitrogen orbitals. 
The highest 4rf-dominated state is unoccupied, consistent 
with the 4d 8 Pd(II) electronic configuration. Adsorbing 
PdP on to Au(lll) perturbs these Pd 4d-like orbitals by 
only a small fraction of an eV, and their occupancies do 
not change. 

To summarize, using the PBE exchange correlation 
functional, we find that PdP lies flat on Au(lll) surface, 
with a binding energy of 0.27 eV. It does not adsorb in 
an edgewise geometry. This conclusion is in qualitative 
agreement with experimental observation of substituted 
porphyrins lying flat on gold surfaces AiSiiiiiS The inter- 
action is dominated by van der Waals forces and is not 
surface site specific. 

D. Mn(II) porphine on Au(lll) — PBE 
predictions 

The adsorption of MnP on Au(lll) is qualitatively dif- 
ferent from that of PdP. First, we consider PBE predic- 
tions using 4 Au layers and T-point sampling. While PBE 
does not predict the correct MnP high spin ground state, 
we will show that it yields qualitatively correct trends 
for MnP adsorbed on Au(lll) in the absence of an elec- 
tric field. Since DFT+U calculations are more expensive 
than PBE, it is also more convenient to use the latter for 
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FIG. 3: (a) Distribution of electron density in a y-cross sec- 
tional plane through the Pd atom. Light (dark) shaded re- 
gions denote densities of at least 0.1 (1.0) electron/A 3 , (b) 
Electronic density of states for palladium porphine (PdP) ad- 
sorbed at the Au(lll) top site. The solid and (inverted) 
dashed curves depict the DOS of the PdP-Au complex and 
of an isolated PdP molecule, respectively, broadened by gaus- 
sians of width 0.1 eV. The sharp lines depict contributions 
from orbitals which have substantial Pd 4d character. The 
Fermi levels (Ef) are shown as dotted lines. 

convergence tests with respect to system size. 

Table |21 shows that the PBE binding energies for MnP 
on Au(lll) are significantly larger than for PdP, and they 
are more site-specific. The top, bridge, and hollow site 
-Ebind differ by up to 0.12 eV. The optimized geometries 
indicate that Mn ion is displaced out of the porphine 
plane towards the gold surface by more than 0.2 A. On 
the other hand, the porphine ring remains ^ 3.5 A from 
the gold surface, demonstrating that the stronger attrac- 
tion arises from the Mn ion interacting with the gold sub- 
strate. The magnitude of the out-of-plane Mn displace- 
ment towards Au(lll) correlates with increasing binding 
energy as the adsorption site varies. When a H2O ligand 
is bound to the Mn ion at the top site, the Mn ion be- 
comes almost coplanar with the porphine ring, and -Ebind 
is slightly reduced, to a value almost identical to that of 
the unligated MnP at the hollow site. In contrast, an iso- 
lated MnPCl is already in the Mn(III) oxidiation state, 
and its adsorption energy is drastically reduced. 

Table 0] shows that 4 layers of gold atoms and T- 



site -Ebind m ri N d p cLau rf p /A 



top 


-0.800 3.77 2.010 -0.224 3.27 


3.49 


bridge 


-0.740 3.70 2.011 -0.216 3.28 


3.50 


hollow 


-0.680 3.68 2.014 -0.205 3.29 


3.50 


top(H 2 0) 


-0.680 3.76 2.015 +0.005 3.55 


3.54 


top(Cr) 


-0.179 3.74 2.028 +0.242 3.74 


3.50 



TABLE 3: Flat adsorption geometry binding energies (eV), 
total spin magnetic moment in the supercell (/ib), and ge- 
ometries of MnP, Mn(H20)P, and MnClP adsorbed at various 
Au(lll) sites. Distances are in A. The PBE exchange corre- 
lation functional, 4 layers of Au atoms, and P-point Brillouin 
zone sampling are applied in these calculations. The symbols 
are described in the caption to Table |5] 



point sampling converge the binding energy to within 
~ 0.09 eV. i?bind does not convergence monotonically 
with system size, and depends on the degeneracy of the 
highest unoccupied orbitals of the Au(lll) slab. Using T- 
point sampling, the same Fermi level degeneracy repeats 
itself whenever 3 layers of gold is added, and, -Ebind are 
almost identical for the 3- and 6-layer gold models. With 
more Appoint sampling in the x-y directions, better con- 
vergence is obtained. We will adopt the 4-layer system 
with T-point sampling as a compromise between accuracy 
and computational convenience. Table 0] emphasizes the 
importance of proper Brillouin zone sampling when treat- 
ing adsorption of transition metal embedded molecules 
on metal surfaces. It suggests that using a finite sized, 
gas phase cluster geometry to represent the gold sub- 
strate, which is always limited to T-point sampling and 
results in an electronic insulator, may be especially prob- 
lematic. 

The spin magnetic moment (to) of the simulation cell is 
also converged to ~ 0.2 /is with 4 layers of Au. We have 
adopted the convention, frequently used in the DFT+U 
literature^*^ that reports the spin magnetic moment m 
as the difference in occupation numbers between up-spin 
and down-spin orbitals, in units of ^b- The magnetic mo- 
ment changes from the PBE prediction of m = 3/ib for 
isolated MnP to M = 3.8/ib for the supercell containing 
the adsorbed molecule. Decomposition of spin densities 
on atomic centers show that the magnetic moment is lo- 
calized on Mn 3d orbitals. Note that, within the frame- 
work of spin-polarized density functional theory, adding 
an electron to the periodically replicated gold substrate 
contributes zero net spin to the entire system. 

In this work, we have adopted isolated ("gas phase") 
Mn(II)P as the reference to compare with MnP adsorbed 
on Au(lll). Although the Mn(II) oxidation state is gen- 
erally less stable in porphyrins than Mn(III), proper- 
ties of stable substituted Mn(II) porphyrins have been 
measured^ and electrochemically induced transitions 
between Mn(II) and Mn(III) are readily realizable. 1,83 
In reality, a gold electrode represents an infinite elec- 
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Q 


1 V/ 1 V/ 1 

1 X 1 X 1 


-0.869 


3.83 


n 

6 


O w O w 1 

Z X Z X 1 


-0.688 


3.19 


6 


o x o x i 


-0.727 3.58 


4 


lxlxl 


-0.800 


3.77 


4 


2x2x1 


-0.783 


3.79 


5 


lxlxl 


-0.693 


3.48 


5 


2x2x1 


-0.786 


3.57 


6 


lxlxl 


-0.868 


3.80 



TABLE 4: Convergence of MnP binding energy (eV) and 
spin magnetic moments (^b) at the top site of Au(lll) as the 
number of Au layers and Brillouin zone sampling vary. 

tron reservoir; when a non-ligated MnP adsorbs on gold 
electrode, the final spin state, charge configuration, and 
molecular geometry will not be affected by the initial 
MnP charge state, although the binding energy measured 
from the Mn(III)P reference state would be different. 

E. Mn(II) porphine on Au(lll) — DFT+U 
predictions 

In this section, we apply the DFT+U technique to ex- 
amine MnP adsorption at the top site of Au(lll). We use 
4 layers of Au atoms with T-point sampling, previously 
shown to yield converged -Ebind from PBE calculations. 
DFT+U consistently predicts a larger -Ebind for the top 
site than the bridge site, just like the PBE exchange cor- 
relation functional. 

Figure Q] depicts the electronic density of states in 
the two spin channels of both isolated MnP and the 
MnP-Au(lll) complex. As already shown in Tabled 
DFT+U (unlike PBE) yields a to = 5^ B (S = 5/2) MnP 
ground state, in agreement with experiments for substi- 
tuted MnP. All majority spin orbitals with substantial 
Mn 3d character are occupied in isolated MnP, while no 
such minority spin orbitals are occupied. 

Adsorption of MnP on Au(lll) leads to striking 
changes in the spectra of orbitals with substantial Mn 3d 
character, which are perturbed by up to 2 eV. The 
highest occupied majority spin Mn-3d orbital in isolated 
MnP is pushed into the conduction band. Overall, the 
MnP-Au(lll) complex exhibits a magnetic moment of 
to = 3.9 hb, similar to PBE results. While both DFT+U 
and PBE predict similar adsorption geometries and bind- 
ing energies (see Tables |3] and their respective Mn 3d- 
like orbitals (not shown for PBE) differ in energies by up 
to 2 eV. 

Figure|SJi illustrates the charge densities in a cross sec- 
tional y plane through the Mn ion and the first layer Au 
atoms. In contrast to PdP (Fig. [3J , there is significant 
overlap of charge densities between Mn and the top site 
Au atom closest to it. On the other hand, the 7r-electrons 




E(eV) 

FIG. 4: Same as Fig. EHb), but for MnP computed using 
the DFT+U exchange correlation functional. The upper and 
lower panels refer to the majority and minority spin channels, 
respectively. The sharp lines are states with substantial Mn 
3d character. 



species 


field 


-Bbind m 


c!n dp dAu 


dp/Au 


MnP 


0.000 


-0.686 3.87 


2.022 -0.247 3.25 


3.50 


MnP 


0.700 


-0.145 4.88 


2.156 -0.540 2.94 


3.48 


MnP* 


0.700 


-0.152 3.22 


2.011 -0.125 3.39 


3.52 


Mn(H 2 0)P 0.000 


-0.737 3.91 


2.020 +0.005 3.55 


3.55 


Mn(Cl)P 


0.000 


-0.106 3.94 2.040 +0.187 3.68 


3.49 



TABLE 5: The binding energies (eV), the total spin mag- 
netic moment in the supercell (^b), and molecular geometries 
of MnP and Mn(H 2 0)P in zero and 0.7 V/A applied field. 
Distances are in A. The DFT+U method is applied except 
for the case marked with the asterisk, which uses the PBE 
functional. The symbols and system size are described in the 
caption to Table |3] 



on carbon atoms of MnP clearly do not interact with 
Au(lll). Comparing the charge density of isolated MnP 
(Fig. [5J;) with that of adsorbed MnP, the top side, out- 
of-porphine-plane electron density in the isolated MnP 
is diverted to the substrate side, where it overlaps with 
gold orbitals. 

To determine the extent of electron transfer between 
MnP and Au(lll), we make a cut through the narrowest 
region of the funnel-shaped electron density distribution, 
and assign charges to the two species accordingly. We 
find there is at most a transfer of 0.2 electron from MnP 



9 




o o au c> 



FIG. 5: Distributions of electron density in a y-cross sectional 
plane through the Mn atom, computed using the DFT+U 
method, (a) No applied field; (b) 0.7 eV applied field; (c) iso- 
lated MnP, zero field; (d) electron density difference between 
(a) and (b). In panels (a), (b), and (c), light (dark) shaded 
regions denote densities of at least 0.1 (1.0) electron/A 3 . In 
(d), the lighter (darker) shades depict regions of electron loss 
(gain) by at least 0.2 electron/A 3 upon applying the field. 
The circles drawn with solid (dashed) lines donote the Mn 
and Au atom positions with (without) applied field. Note 
that these atoms and their electron clouds move in opposite 
directions. 



to the substrate, despite the drastic changes in the DOS 
and the spin multiplicity from their gas phase values. The 
binding between MnP and Au(lll) is thus best described 
as covalent or metallic. 



F. Applying an electric potential 

Given that the MnP-Au(lll) interaction is associated 
with partial electron transfer, we expect that changing 
the electric potential on the gold substrate will strongly 
influence the adsorption behavior. An electric field 
which favors transfer of electrons from MnP to Au can 
strengthen the MnP-Au interaction. An electric field 
which favors electron transfer back into MnP 3d orbitals 
should cause the Mn(III)-like behavior to revert back to 
Mn(II). 

Applying electric fields within plane wave density func- 



tional theory calculations has been formulated in the lit- 
erature and implemented into the VASP code^ We ap- 
ply a field of 0.7 V/A to isolated MnP, a 4-layer gold 
slab, and the adsorbed MnP on 4-layers of Au atoms. 
£"bind is obtained as the energy difference between the 
adsorbed and isolated systems. T-point sampling is used 
throughout. 

Table[S]shows that a 0.7 V/A electric field strongly re- 
duces the binding energy and changes the magnetic mo- 
ment to almost its isolated MnP value: m — 4.89^tB ac- 
cording to DFT+U predictions, compared with m = 5/ib 
in isolated MnP. The Mn-N distance, d^, also increases 
from the zero field value of 2.02 A to 2.15A, which is even 
larger than the value of = 2. 090 A predicted in iso- 
lated MnP. The PBE predictions for molecular geometry 
and spin states are substantially different; this reflects 
the inability of PBE to predict the correct spin state and 
molecular geometry for isolated MnP. 

The larger Mn-N distance in adsorbed MnP in an elec- 
tric field compared to isolated MnP is due to the sig- 
nificant Mn ion displacement out of the porphine plane 
towards the gold substrate, and not to a change in the 
porphine conformation, which remains essentially planar. 
This is demonstrated by the similar core sizes for isolated 
MnP (2.090 A) and for MnP adsorbed on Au(lll) with 
an electric field (2.087 A) (the core size is defined as the 
radius of a cylinder which can fit through the porphine 
hole) . NSD analysis (supporting information) also shows 
that the A lg (in- plane) deformation of MnP on Au(lll) 
in the absence of an applied field is consistent with a 
Mn(III)-like species, with a Ai g deformation of +0.08 A 
compared to +0.11 A predicted for isolated MnClP. As 
expected, with an applied electric field of 0.7 V/A, the 
first order Ai g deformation of +0.31 A is similar to that 
predicted for isolated MnP (+0.37 A). 

Figure [5] illustrates the effect of the applied field on 
the charge density. The field evidently repopulates the 
Mn 3d x 2_ y 2-\ike orbital and restores electron density to 
the top side of the MnP molecule. There still appears 
substantial electron cloud overlap between MnP and the 
gold surface despite the small £bind- Figure [SJi depicts 
the induced charge differential. Note that the field causes 
both the surface top site Au atom and the Mn ion to move 
towards the gold substrate, whereas the electron densities 
on both atoms are displaced in the opposite direction. 

Figure HJ1 plots the electronic DOS predicted by 
DFT+U in the applied field. Compared with the zero 
field case (Fig.^J), the occupied states are shifted to lower 
energies relative to the Fermi level. In particular, the 
highest lying orbital with substantial Mn 3d character is 
brought below the Fermi level in the MnP-Au(lll) com- 
plex, and is occupied by electrons. This is consistent with 
the field induced change in the total magnetic moment. 
Nevertheless, if we demarcate the charge density at the 
narrowest region of the electron cloud between Au and 
Mn, once again we find a small, ~ 0.1 electron difference 
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FIG. 6: Same as Fig.|U for MnP computed using the DFT+U 
method with an applied electric field of 0.7 eV/A. 



between the total charge on MnP and Au with or with- 
out the electric field. Thus, while the charge and spin 
state of the Mn ion differ, the locus of the total electron 
density is not strongly affected. 

The 4-layer bare gold slab is stable in the 0.7 V/A 
field, and the field-induced change in its energy is pro- 
portional to the square of the field magnitude. We find 
that smaller, 0.525 V/A and 0.35 VA electric fields lead 
to similar behavior described for MnP adsorbed on gold 
in a 0.7V/A field, although the Mn-N bonds becomes 
shorter (2.109 and 2.102 A, respectively). Due to the 
technical difficulty of converging a DFT+U calculation 
in a system with spin polarization, a zero band gap, and 
an electric field, we have not located the minimal field 
needed to trigger this switching. This will be investi- 
gated in a future work. These calculations are performed 
in vacuum, and we have not attempted to relate the field 
strength to an electrostatic potential. In electrochemical 
cells, where the molecules and electrodes are immersed 
in electrolytes, the voltage difference required to reduce 
Mn(III) porphyrins to Mn(II) is much smaller: about 
0.4 Volt in water^ and on the order of 0.2-0.3 Volt in 
organic solvents i 



IV. CONCLUSIONS 

In this work, we use a combination of density func- 
tional theory (DFT) in the generalized gradient approx- 
imation (GGA) and the DFT+U technique to study the 



adsorption of transition metal porphine molecules on 
atomistically flat Au(lll) substrates. The PBE exchange 
correlation functional is adequate for treating palladium 
porphine (PdP). We find that PdP preferentially adsorbs 
flat on Au(lll) surfaces, in agreement with experimental 
observations of substituted porphyrins lying flat on gold 
surfaces AiMiii^ The binding energy is 0.27 eV, irrespec- 
tive of the adsorption site. There is no charge transfer 
or covalent bonding between PdP and Au, and the inter- 
action is predominantly dispersive. In contrast to a pre- 
vious DFT calculation^ we find that a self-assembled 
tube-like geometry with the PdP adsorbing edgewise on 
Au(lll) is not energetically favorable. Using multiple 
DFT codes and types of pseudopotentials, we find that 
the LDA binding energy for this edgewise adsorption ge- 
ometry is 0.25 eV, not ~ 10 eV as previously predicted^ 
The PBE exchange correlation functional predicts little 
or no binding. Furthermore, an accurate binding energy 
for this geometry requires more extensive Brillouin zone 
sampling than has been attempted in the literature. 

The behavior of manganese porphine (MnP) is qualita- 
tively different. The PBE exchange correlation functional 
does not reproduce the correct ground state (high spin 
5 = 5/2 state) or the experimental Mn-nitrogen distance. 
However, PBE can be augmented with the DFT+U tech- 
nique, which can be parameterized to yield the correct 
ordering of spin multiplicity. The DFT+U ground state 
geometry exhibits a Mn-N distance in good agreement 
with experiments and B3LYP predictions. 

MnP selectively binds to the top site of Au(lll). When 
we use Mn(II)P as the reference state, the binding energy 
of MnP on Au(lll) is significantly larger than for PdP. 
The adsorbed MnP exhibits partial charge transfer from 
orbitals with substantial Mn 3p character to Au(lll) 
states. The overall magnetic moment also changes from 
m = 5/ib (S = 5/2) to m « 4/ib (5 = 2). In other words, 
Mn takes on a +3 oxidiation state character. 

This MnP-Au(lll) interaction can be partially re- 
versed by applying an electric field with the appropri- 
ate polarity. A sufficiently large field causes the Mn ion 
to revert back to a Mn(II)-like spin state, resulting in a 
significant, up to 0.13 A, increase in the Mn-N distance 
within the porphine ring. This suggests that appropri- 
ately substituted Mn porphyrins deposited on gold elec- 
trodes may be useful for demonstrating electric field trig- 
gered conformational changes that are potentially per- 
tinent to harvesting nanomechanical work and selective 
ligand binding^ 

In summary, we have successfully applied the DFT+U 
technique to transition metal ions in a condensed phase 
environment with fully periodic boundary conditions. 
With this approach, the energy difference between the 
spin states can be parameterized using not just B3LYP, 
which gives qualitatively correct spin orderings for many 
manganese complexes, but experimental data, quantum 
chemistry methods, and quantum Monte Carlo calcula- 
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tions as well. As such, this versatile technique can poten- 
tially be used for general, accurate, all DFT-based treat- 
ment of transition metal ions in aqueous systems, water- 
material interfaces, and biological environments M£k 
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